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The effective chiral interaction between molecules arising from long-range quantum interactions 
between fluctuating charge moments is analyzed in terms of a simple model of chiral molecules. This 
model is based on the approximations that a) the dominant excited states of a molecule form a band 
whose width is small compared to the average energy of excitation above the ground state and b) 
biaxial orientational correlation between adjacent molecules can be neglected. Previous treatments 
of quantum chiral interactions have been based on a multipole expansion of the effective interaction 
energy within second-order perturbation theory. We consider a system consisting of elongated 
molecules and, although we invoke the expansion in terms of coordinates transverse to the long axis 
of constituent molecules, we treat the longitudinal coordinate exactly. Such an approximation is 
plausible for molecules in real liquid crystals. The macroscopic cholesteric wave vector Q (Q = 2n/P, 
where P is the pitch) is obtained via Q = h/K2, where K2 is the Frank elastic constant for twist and 
h is the torque field which we calculate from the effective chiral interaction kijSli x a./ • Ru, where 
the unit vector aj specifies the orientation of molecule I and R/j is the displacement of molecule I 
relative to molecule J. We identify two distinct physical limits depending on whether one or both 
of the interacting molecules are excited in the virtual state. When both molecules are excited, we 
regain the Rjj dependence of kij on intermolecular separation found previously by van der Meer et 
al. The two-molecule, unlike the one-molecule term, can be interpreted in terms of a superposition 
of pairwise interactions between individual atoms (or local chiral centers) on the two molecules. 
Contributions to kij when one molecule is excited in the virtual state are of order Rjj for helical 
molecules which are assumed not to have a global dipole moment, but whose atoms posses a dipole 
moment. It is shown that for a helical molecule Q can have either the same or the opposite sign as 
the chiral pitch of an individual molecule, depending on the details of the anisotropy of the atomic 
polarizability. The one-molecule mechanism can become important when the local atomic dipoles 
become sizable, although biaxial correlations (ignored here) should then be taken into account. Our 
results suggest how the architecture of molecular dipole moments might be adjusted to significantly 
influence the macroscopic pitch. 
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I. INTRODUCTION. 



In the cholesteric liquid crystalline phase,™ 2 ] anisotropic mesogens align on average along a local unit director n(r) 
that rotates in a helical fashion about a uniform pitch axis. The pitch P of this helix ranges from a few tenths of a 
micron jtftten or more microns. In fact, solutions of the viruses FD and TMV, as well as DNA, have even much larger 
pitches.tffl Because the pitch is usually large compared to the intermolecular separation, these systems are locally 
essentially indistinguishable from nematics and consequently they are often referred to aSptjiiiral nematics (CN's). 
The pitch wavenumber Q = 2n/P can even pass through zero as a function of temperature. EiEl The helical structure 
of a cholesteric phase must result from the molecular chirality of some or all of its constituent mesogens. Achiral 
mesogens form an achiral nematic rather than a chiral nematic phase. Phenomenologically, the explanation of the 
twist of the cholesteric phase is straightforward: chiral mesogens must lead to a chiral term hn • V x n in the long- 
wavelength free energy density that favors twist. This tendency to twist is resisted by a twist elastic energy density 
ii^2(n -V x n) 2 , where K 2 is the Frank elastic constant for twist. If the pitch axis coincides with the z-direction, 
then in the equilibrium configuration one has 

n(r) = (cos Qz, sinQz, 0) , (1) 

with 



Q = 2ir/P = h/K 2 . (2) 

(Our definition of Q is such that positive Q corresponds to right-handed macroscopic chiralityfi) The magnitude of K 2 , 
which has units of energy per unit length, is estimated with good accuracy by dimensional analysis. The characteristic 
energy is of order the thermal energy, kgT ~ fceT/v/ where kg is the Boltzman constant, T is the temperature, and 
Tjvj is the isotropic-to-nematic transition temperature. The characteristic length is a molecular length L, so that 
K2 ~ ksT/L. A similar dimensional analysis for the torque field h, which has units of energy/ (length 2 ), would predict 
h ~ kgT/L 2 and P ~ L. This is a far tighter pitch than is observed in any cholesteric. This reasoning indicates that 
an explanation of the magnitude of h requires considering a detailed model of the cholesteric. The chiral structure of 
cholesterics also raises some technological issues. It would be very desirable to be able to "engineer" molecules that 
have specific values of h and thus P or more generally that have a specific temperature dependence for h. To realize 
this goal, it is necessary to understand how variations in molecular architecture and electronic structure influence h. 
As a first step in dealing with these issues of fundamental and applied science, this paper will address some aspects 
of the calculation of h from a molecular model. 

In fact, the calculation of h is highly nontrivial. It involves the rather complex interactions between mesogens 
and the orientational correlations they induce. If there are chiral mesogens, there are chiral interactions, and h 
is nonzero; otherwise, h is zero. One typically identifies three types of interactions between moleculesfa (1) long- 
range attractive dispersion (Van der Waals) interactions, (2) short-range repulsive interactions, whose origin is the 
Pauli principle, and (3) direct Coulomb interactions, which take the form of dipolar, quadrupolar, etc. interactions 
between electrically neutral mesogens. The latter interactions are of secondary importance in many chiral and achiral 
liquid crystals and will be ignored here. Initially Straleyfl proposed that the macroscopic chirality of CN's could be 
understood qualitatively n i]x-|jeijms of the packing of screws. These short-range repulsive forces were modeled as hard- 
core or steric potentials,ETt£rE3 reflecting molecular shape, that contribute to the entropy but not the internal energy. 
For spherical atoms, the repulsive and dispersion forces can be combined in a single effective potential such as the 
Lennard- Jones 6—12 central-force potential. More generally, interactions between achiral molecules can be modeled 
as sums over central-force effective potentials between pairs of atoms or mass points on different molecules. til There 
are chiraL-wersions of both dispersion and short-range xepulsive forces. Chiral dispersion forces were first analyzed by 
GoossensEa and later more systematically by others .fl>E3EJ They found that the dominant chiral interaction between 
chiral mesogens, calculated in the limit of center-of-mass separation R much larger than any molecular dimension L 
was proportional to R~ 7 and to the product of dipolar and quadrupolar molecular matrix elements. Various^pmewhat 
ad hoc chiral intermolecular interactions, some based, on implementing models equivalent to threaded rods,Ejll3 others 
on surface-nematic interactions of chiral dopants,tj have beea, introduced mostly as input to simulations of chiral 
systems. Models for flexible mesogens have also been treated.c2l I— . 

A chiral molecule is one that cannot be rotated into coincidence with its mirror imageEd Chiral molecules cannot 
be uniaxial: at minimum, their description requires an orthonormal triad of vectors rather than a single vector. A 
microscopic description of chiral interactions involves the complete orthonormal triad of axes emblazoned on each of 
the two interacting molecules. However, as we have mentioned, apart from very small corrections arising from slow 
local twist, the cholesteric phase is locally uniaxial. It is, therefore, natural to seek effective chiral interactions between 
effectively uniaxial molecules. If a molecule J of arbitrary shape is spun about some axis aj, it becomes on average 
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uniaxial with respect to this axis. Thus, general pair interactions between molecules / and J in a chiral nematic can 
be reduced to uniaxial pair interactions by averaging over independent rotations of each member of the pair about the 
local nematic director. The resulting potential is only approximate in that it ignores orientational correlations between 
molecules in the plane perpendicular to aj and slj. In practice, one usually averages over independent rotations of 
each molecule about its body axis a, rather than the more correct average over rotations about the local nematic 
director. We mention that it is known that thej-chital part of central-force potentials (such as hard-body interactions) 
vanishes when such correlations are neglected PW^ However, dispersion contributions to the chiral interaction do not 
require nonzero orientatipaal correlations between molecules. As discussed in Appendix ^ and as has been found by 
several previous authors ) 13 H 14 I the long-range dispersion interaction survives this independent rotation procedure to 
produce an effective chiral potential between effectively uniaxial mesogens of the form 

#7j(aj,aj,R/j) = (a/ x a 7 • Ru)ku, (3) 

where R/,/ is the displacement of the center of molecule / relative to the center of molecule J and only terms in kij 
which are odd in both slj and slj are retained. The effective interaction of Eq. (||) arises between two chiral molecules 
as well as between a chiral molecule and an achiral one. . .. ■ 

Our derivation of the effective chiral interaction differs from previous onescJO in two important respects. First, 
previous calculations of this interaction are based on a multipole expansion in the variable Ti/Rjj, where is the 
coordinate of the zth charge of molecule I relative to the center of molecule I. Strictly speaking, the multipole 
expansion only applies when R/j = R is large compared to any dimension of the molecules. This expansion does 
not apply to pairs of molecules whose separation is less than their length but greater than their width. We develop a 
modified multipole expansion in which coordinates transverse to long molecular axes are treated as small parameters. 
Second, the results of previous calculations are expressed in terms of electric dipole and quadrupole matrix elements 
of the entire molecule. But in .a-Jong molecule, typical of those comprising liquid crystals, we expect the electronic 
states to be strongly localizedo Accordingly, it seems more useful to express results in terms of matrix elements 
within atoms or local complexes. In so doing, it is natural to assume the relevant excited state can be reached from 
the ground state by matrix elements of the dipole moment operator. Then, the quadrupole moment operator is 
easily related to the dipole moment operator, with the result that the only matrix elements appearing in the present 
paper are those of the dipole moment operator between local atomic states. In common with previous treatments, 
we will neglect the effects of biaxial correlations between interacting molecules. Accordingly, we. will evaluate kij by 
averaging each molecule independently over spinning about its long axis. In a separate paperc3 we will discuss how 
the chiral interaction between helical molecules depends on the angles describing rotation about their longest body 
axis. 

We may summarize briefly the results of this program. Although we do not expand in powers of the longitudinal 
coordinates of the charges in each molecule, our results are formally not very different from the previous oneslldca 
However, by expressing the results in terms of matrix elements of localized atomic orbitals, we identify two distinct 
physical mechanisms. The first is the dipole-quadrupole interaction previously identified. The second is one involving a 
three-body interaction between two local atomic dipole moments on one molecule and a local anisotropic polarizability 
of the second molecule. This second interaction, formally present in previous work, can dominate the first one in 
certain situations. Furthermore, our approach allows us to discuss how these interactions depend on the length of 
the molecule. For a helical molecule we find that the contribution to kjj due to the first mechanism is proportional 
to L 2 /R s for L ^ R and to L/R 7 for L > R. Results for the three-body interaction are more complicated [see Eq. 
(|65|), below]. In both cases, the magnitude of the pitch arising from these interactions in a concentrated system of 
helical molecules with polarizability corresponding to a dielectric constant of about 1.3 would be 10 microns. This 
is a larger pitch than one observes for most concentrated cholesterics. It is possible that this discrepancy is due to 
some of the, simplifying assumptions in our calculations most probably our disregard of biaxial correlations between 
moleculesnJ Alternatively, it is possible that the pitch of most cholesterics is determined by steric rather than by 
quantum interactions. Elsewhere we will apply the approach of the present paper to obtain the quantum contributions 
to K 2 !B 

Briefly, this paper is organized as follows. In Sec. II we give an overview of the calculation of the torque field h 
from which the macroscopic chiral pitch can be determined. In Sec. Ill we derive a rather general expression for 
the strength of the chiral interaction between molecules in terms of matrix elements of the dipole moment operators 
of atoms between the ground state and excited states localized on atoms. The contribution to this effective chiral 
interaction from virtual states in which both molecules are excited is treated in Sec. IV and that in which only one 
molecule is excited is treated in Sec. V. Numerical estimates of the pitch for a system of helical molecules are given 
which show that the quantum mechanism when both molecules are excited is unlikely to explain the observed pitch 
of most cholesterics. Our results and conclusions arc summarized in Sec. VI. 
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II. OVERVIEW OF CALCULATION 



In this section we give an overview and a summary of the results of the calculation, given in the next section, of the 
chiral interaction between molecules. Recentfyo a systematic formulation was given that expresses the macroscopic 
pitch of a CN in terms of microscopic interactions between molecules. Such a formulation is required in cases where 
it is either necessary or desirable to include orientational correlations between interacting molecules. It was shown 
that for central-force interactions between atoms on different molecules, a nonzero effective chiral interaction between 
molecules could only be obtained when orientational correlations, specifically biaxial correlations, between molecules 
were taken into account. In contrast, in the present paper we will see that the quantum interactions between molecules 
are not of this type. Thus, in the present context, it is permissible to use a simpler and more traditional approach 
in which each molecule is characterized by the orientation of its long axis, specified by the unit vector a. We will 
then evaluate the chiral interaction energy between molecules I and J written in Eq. (|J). This interaction energy 
is evaluated within what we will call the "uniaxial" approximation in which we independently average over the 
orientations of the two molecules when a is specified for each molecule. In order to calculate the macroscopic chiral 
pitch it is only necessary to evaluate the chiral part of this interaction, i. e. the part of the form written in Eq. 
(|J) . This interaction energy can then be added to whatever phenomenological interaction one is using to describe the 
nematic phase which would result in the absence of chiral interactions. 

To make contact with a continuum theory, one introduces a local order parameter tensor via 

Q Q/3 (r) = a a (r)ap(r) - ±<5 Q/3 . (4) 

When thermally averaged, this tensor becomes the usual de Gennes-Maier-Saupe order parameter. In the long 
wavelength limit, the chiral interaction Eu leads to a continuum interaction of the form 

Eh& = i J dre af ^Q aS (r)V f3 Q lS {r) , (5) 

where Greek indices label Cartesian components, 6 a p is the Kronecker delta, e Q ^ 7 is the antisymmetric tensor, the 
repeated-index summation convention is understood, and the constant 7' is the macroscopic analog of kij. If we 
express Q in terms of the director n(r) as 

(QapW = S (n a np - \8 a p) , (6) 

where ( )t indicates a thermal average, then the above chiral interaction leads to th&.familiar Frank free energy in 
the presence of macroscopic chiral twist (but neglecting splay and bend distortions) asB 

F =\K 2 J dr[n(r) • V x n(r)] 2 + hj drn(r) • V x n(r) , (7) 

where, within mean-field theory, h — j'S 2 . The main goal of the present paper is to obtain h from a microscopic 
model. In particular we will obtain h from an evaluation of kij in Eq. (|J). As we have seen in Eqs. (Q) and a 
determination of h leads immediately to the determination of the macroscopic chiral wave vector Q. In addition, the 
chiral properties of an isotropic liquid consisting of chiral molecules, such as the rotary power, are also related to kjj. 

The interaction between molecules we are going to study is the generalization, for chiral molecules, of the attractive 
1/i? 6 term in the van der Waals potential between neutral spherical atoms. This calculation is based on a quantum 
mechanical treatment of the total Coulomb interaction TLu between charges on the two interacting molecules, I and 
J: 

where i G I indicates that the sum is over all charges qi, both electronic and nuclear, in molecule / and Rij is the 
displacement of qi relative to qj . In this calculation we nes^ect any biaxial correlations between the orientations of the 
two molecules. Within this assumption it has been shownE3 that central force interactions [like those of Eq. (|[) when 
taken in first order perturbation theory] can not lead to any chiral interactions. Therefore, we consider the effect of 
H.U within second-order perturbation theory. The effective interaction between molecules I and J is then 

Eu = Eufa^-Rij) = (0\HujHu\0) , (9) 
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where £ — Eq — Ho, with Eq the energy of the ground state |0) of TCq the unperturbed Hamiltonian describing 
noninteracting molecules, and V is a projection operator that excludes the ground state. Here we have indicated that 
Ejj depends on uj (and ujj) which denotes the triad of Euler angles (see Fig. |]) needed to specify the orientation 
of the Jth (Jth) molecule and on R/,/. Both R/j and r^, the position of the ith charge in molecule I relative to the 
center of the molecule, r^, may be expressed with respect to coordinate axes, e M , fixed in space, as shown in Fig. ||: 

R/ = i?/ Al e Al , Ti = r^e^ . (10) 

We now must average this interaction energy over the orientations of the two molecules, when the long axes of 
molecules I and J are fixed to lie along the unit vectors a/ and aj, respectively, and correlations between the 



orientations of the two molecules are neglected. To carry out this average we introduce axes specified by unit vectors 
e' Ifl emblazoned on the 7th molecule, as shown in Fig. [y, so that 



TV i 

I 

n = ^V e /M ■ ( n ) 

Previously van der Meer et al.lll carried out this averaging within the multipole expansion. However, since we wish to 
treat long molecules of the type usually constituting liquid crystals, we do not make the usual multipole expansion, 
but rather expand only in terms of transverse coordinates of the molecule. Thus we set 

r t = z'pi + p t , (12) 

where p i ■ aj = 0. [Throughout, atom i (j) is assumed to be in molecule I (J)]. Thus z\ and p i are the coordinates 
of the ith charge of molecule I relative to the center of the molecule, respectively, longitudinal and transverse to the 
long axis of the molecule aligned along a/ = e' z . Now we expand Ejj in powers of p and perform the orientational 
average over powers of p (indicated by [ ] av ) is done using, e. g., 

[(pi)a(Pi')p\*v = U x i x i' + y'iVi'){ 5 OL,0 - a a a ) + |e a/ 3 7 a 7 (x-y-, - y>-,) 

= l^i'^OW - ciaap) + ^e^^r^ea^a^ , (13) 

where /j, and v run over only transverse (x, y) coordinates and a 7 = (a/) 7 . Thereby we find that 

[-B/j]av = (a/ x a, 7 • R/j)K/j(a/, aj, R/j) + . . . . (14) 

Here we have written the term responsible for the chiral interaction between molecules / and J and have discarded 
the nonchiral terms (represented by . . .). 

The expressions for kjj in its most general form are not very enlightening, although they do display the appropriate 
symmetry to vanish for molecules which are not chiral. To gain some insight into the meaning of these results we have 
had recourse to a model of the excited states, which appear in second-order perturbation theory. Our first assumption 
is that the important excited states consist of dipolar fluctuations from the ground state. In other words these states 
are taken to be the three atomic p states on atom i. The second assumption is that 5, the width in energy of the 
band of excited states obtained by allowing these excitations to occur on any atom is small compared to their energy E 
relative to the ground, state. This assumption allows us to take the virtual intermediate states to be strictly localized 
to individual atomsBJ Nonlocal effects give rise to corrections of relative order, 8/E. Under these assumptions, our 
results may be summarized as follows. Contributions to kjj can be classified into two types, depending on whether 

one or both molecules in the intermediate state are in an excited state. These are denoted k\ j and nfj , and will be 
referred to as "one-molecule" and "two-molecule" terms, respectively. Our results are conveniently written in terms 
of the definition Kj 1 ) = ^[ft/j + ]a ; where, for any function / of &j and aj, 

[/] a = i[/(a / ,a J )-/(-a J ,a J )-/(a / ,-a, 7 ) + /(-a / ,-a J )] . (15) 



Then 



where 



r.(2) 



S t , = [a, • a, - 2( a/ • D«)(aj • T),J if, I),; (17) 
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Mij = 3e 4 



E 



% <0| Az'j \ H ) (nj\Ax'AO) - JBj (0| A^-X^l A^|0> 



2(0|Az^) 2 - (0|Ax' ( |^) 2 - (OIA^I^ 



(18) 



where the sums over i and j now run only over electrons, r'j is the expectation value of r^- in the ground state, i. e. 
it is the center of the atom associated with charge j, Ar'j = r'j = r'j and 



R 



- z^a/ - z 1 , aj 



(19) 



Also \fij) denotes the state when all atoms are in their ground state except for atom j, which is in the excited p state 
labeled fi, which has energy E^{j) relative to the ground state and Ey^{i,j) = E„{i) + E^(j). Here | Ar^ |0) is 
nonzero only when j refers to an electronic charge. We also find that 



k {1) 
K u 



ii'j 
„2 ^ "" > 



[PixPi'y Pi'xPiy] 



(ggj ■ g£) 



x e 2 \ ^OrMs^lAz'lO) 2 - ( H \Ax'M 2 - (/^A^O) 2 ] 



(20) 



where the sums over i and i' are over atoms and p' Q is the a component of the dipole moment vector in the ground 
state evaluated in the molecule-fixed coordinate system. (For this calculation a local dipole moment was assumed, 
but this does not necessarily imply the existence of a dipole moment of the molecule as a whole.) Note that the above 
expressions, since they have already been averaged over rotations about the long axis, are invariant with respect to 
rotation of each molecule about its long axes parallel to e' z . 

Our result for k^J is closely related to that previously obtained by van der Meer et al@ and by Kats.0 To obtain 
a form close to that obtained by Kats, we write 



<j'(j;u) + iQ + )^'{i\-uj + i0 + )dui , 



(21) 



where 



A* 



(22) 



</(*; u) = e 2 (2(0\z<W 3 f - {0\x^f - (0|y>,-) 2 ) /[u - E u {{)] 

= a' zz (i;uj) - \a xx {i\uj) - \a' yy {i,u) , (23) 

where a' is the [i-v component of the polarizabihty tensor with respect to the molecular frame. Here 7'(w) is the 
anisotropy of the polarizabihty and <j'(u>) is the higher order quadrupole-dipole response function (which Kats calls 
the gyrotropy), both taken in the molecular frame, as indicated by prime superscripts. Here these quantities are given 
by a sum over the corresponding properties for the individual atoms. We assume that the relevant excited states are 
localized p states, in which case the gyrotropy can be related to the polarizabihty: 

a'(j;w) = e 2 ^ (^.(0|4|^)(^-|4|0) -^<0|4| Mi )<^-|^|0)) /[w - E„(j)} 

= - \x' 3 a yz {j- : uj) . (24) 

In addition the factor Sij depends on the z'-component of the position of the ith atom, whereas in the bare multipole 
expansion used by van der Meer and Kats, only R/j appears. Because we do not include the z' coordinate within the 
multipole approximation, we can treat long molecules in an appropriate way, as is reflected in the sum over atoms 
of Sij. One sees that the chirality of the molecule is incorporated in a', which vanishes if the molecule has a mirror 
plane.0 In K^j, the chirality of the molecule / is incorporated in terms like 



G 



T ' = ^[PtxPt'y ~ Pi'xP'iylWi ■ (25) 

ii' 

In the case of classical interactions, it was not possible to construct a third rank tensor of the mass moments which 
was zero for achiral molecules and nonzero for chiral molecules,E3 because such a mass moment tensor was symmetric 
under interchange of any pair of its three indices. Here, however, one sees that a' and r are x,y,z elements of tensors 
that are not symmetric in all indices and that, therefore, can be used as an indicator of chirality. 

It is interesting to evaluate these expressions for some specific geometry of a chiral molecule. For this purpose we 
treat in some detail a helical molecule, patterned after DNA. Then we introduce local atomic coordinates whose axes 
coincide with the axes defined by the local excited p states. We assume that these axes, shown in Fig. ||, are identical 
to those of the tangent, the normal, and the binormal unit vectors, which we call e'^, with (i respectively, z, x, and y, 
so that we can relate the a'^ v to its components a'^ ^ in the local atomic frame as 



>-<>- 



(j) = 2e 2 J2 E^r^^W^'L ' OK> ' e *) = («& ' ' ^K„U) ■ (26) 



One should note the following general points in connection with our results. Firstly, the result in Eq. (|T^) shows 

(2) 

that Kjj can be viewed as arising from a superposition of interactions between local-centers of chirality on one molecule 
with centers of anisotropic polarizability on another molecule. As is well knownjlia this result implies that chirality 
can be induced by the interaction between a chiral molecule and an achiral one that has a local center of anisotropic 



polarizability. In contrast, the result in Eq. (20) is a three-body interaction between two local dipoles on one molecule 
(combined with resulting chiral strength r') with a local anisotropic polarizability of the second molecule. Finally, we 
mention that it is interesting to generalize these results to a flexible polymer the orientation of whose backbone may 
vary appreciably over its length. 

III. CHIRALITY FROM INTERMOLECULAR INTERACTIONS 

We now turn to the calculation of kjj. For this purpose we give a brief discussion of how the average over 
orientations is to be done. In general, the orientation-pf the 7th molecule is specified by the three Euler angles or, 
/?/, and 7/, for which we adopt the definition of Rose,c3 as is illustrated in Fig. [l| In particular, ai and fii are taken 
to specify the orientation of the long axis of the molecule. So we write 

slj = sin (3 j cos aje x + sin [3j sin aje y + cos (3je z . (27) 

Within the spirit of mean-field theory we should average the interaction energy between molecules I and J over the 
single-molecule orientational distribution function appropriate to a nematic, which locally is a good description of 
the CN. For molecule I this average should be taken subject to its long axis being specified by the fixed value of aj. 
The single molecule orientation distribution function p(uji) must be independent of aj and also should be invariant 
under a — * —a. Specifically, correlations between j3j and 7/ are allowed^J as is discussed in Appendix 0. For any 
function of molecular orientation /(<£>) this average is 

[/(&)]av = f p(&i)f{&i)d'yidai/ /"p(cjj)d7j . (28) 

However, when the molecule is not strongly biaxial, or when the molecule is. jaerfectly aligned along the nematic 
direction, the assumption that p(u>) is independent of 7/, as is usually done,Ej'E-3 is sufficient. This approximation, 
which we call the uniaxial approximation, will be used in this paper. In addition, to preserve invariance under 
a/ — > —a/, we will also project out of the calculation terms in kjj that are even in a/ and in a,/. This step can be 
done at the end of the calculation using Eq. rUlg) ■ 

Our calculation of kjj, as previous ones,EJE3 is analogous to that of the well-known R~ e interactions between 
widely separated neutral atoms. Quantum fluctuations involving dipole moments in excited states are treated within 
second-order perturbation theory. Short-range quantum repulsion is often treated in an ad hoc fashion via a classical 
central-force interaction between atoms but this effect will not be discussed here. We take the interaction Hamiltonian, 
TLu, for molecules I and J to arise from the Coulomb interaction between the ith charge on molecule /, denoted q^, 
and its counterpart on molecule J. Thus we write 

"'^| i.,:t,r (29) 
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We use Eq. ( [12 ) to write 



2 ^ 1 2 

ij ij 



-1/2 



(30) 



where p^ = p i — p and 



Note that is evaluated for p i = Pj = 0. 

We now expand with respect to transverse coordinates to obtain 
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35 [pjj ■ T> t 
8 ^ 



(31) 



(32) 



Note that this expansion is valid if the charge distributions of the two molecules do not overlap one another. Strictly 
speaking, the validity of our treatment requires satisfying this condition for all configurations with nonnegligible 
weight in the partition function. 

We now consider an evaluation of the interaction energy between two molecules treating TLu via perturbation 
theory. The first term is the ground-state expectation value of the Coulomb interaction between atoms on different 
molecules. If we neglect biaxial correlations between orientations of adjacent molecules and we simply average this 
interaction over the uncorrelated rotations of the two molecules, subject to their long axes being fixed, then we find 
the resulting interaction to be completely achiralE3 Accordingly, to obtain an effective chiral interaction from "ftij 
when such biaxial correlations are neglected, it is necessary to evaluate the energy of interaction within second-order 
perturbation theory, whereby 



E 



ij 



' I (ft Ij) n i, n. j ;0,0 |" 

nj,n j 



(33) 



where the sums are over states |nj) (\nj)) of molecule I (J) and the prime indicates exclusion of the term when 
both molecules are in their ground state. Here E ninj is the energy (relative to the ground state) of the state when 
molecules I and J are in states |nj) and \nj), respectively. 

The obvious step of substituting the expansion of Eq. ([52]) into Eq. ([53]) leads to rather complicated algebra. We 
now classify terms according to their order in Rjj . Since we do not assume the length L of the molecules to be much 
less than the separation Rjj between molecules, in counting powers of Ru we consider L/Rjj ~ 1. As we shall 
see, Eh ~ RJj,c3 where p = 7 for two-molecule terms and p = 6 for one-molecule terms. Accordingly, we drop all 
contributions which are of order RJf with p > 7. Also, we drop contributions which are proportional to odd powers 
of p, since these will vanish when we average over rotation about the long axis of the molecules. Thereby we obtain 
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When one averages over independent rotations of the two molecules about their long axes, using Eq. (|I3|), one sees 
that the first two lines of Eq. (5^) do not lead to a chiral interaction. 

We imagine the virtual statesjVj} in Eq. ( |35| ) to be a linear combination of excited atomic p states. Accordingly, all 
matrix elements can be chosen to be real. Also, in this model we take no explicit account of exchange and correlation 
effects beyond what is included in self-consistent atomic orbitals. Thus, it is permissible to label electrons according 
to their atomic location. Then, for the matrix element of an arbitrary function, /, of we can write 

(m\f(Ti)\0) = V' a /(r) | r<=F (m\Ar' ia \0) + O ((r^Ar^Ar^O)) , (36) 

— Fj. To leading order in l/Ru we have 

.„ , / = -(n i |A4|0)(a J -D y )Sr 3 , 



where Ar^ 
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niUjlAzlAz, \0i0j) [3(ai • D tf )(aj • D y ) 



a 7 • a.jD 2 lj \D l 5 



yj y ' 

(37) 

where |ra,rij) is the state (whose energy relative to the ground state is E nin .) in which atom i is in excited state |rij), 
atom j is in state and all other atoms are in their ground state and was defined in Eq. @. Thus [D^ 1 ] 

is of order at least Rjj for single-molecule terms and of order at least RJ 3 for two-molecule terms. This argument 
shows that the last line of Eq. (|35|) does not contribute at leading order in 1 / Rjj and that we have only to deal with 
lines 3, 4, and 5 of this equation. 

We now carry out the average over the orientations of molecules I and J subject to their long axes being fixed to 
be, respectively, along a/ and aj, using Eq. ([l3]). and we keep only chiral terms of the form written in Eq. (||). This 
procedure is algebraically extremely complicated. However, the fact that the relevant excited states are undoubtedly 
strongly localized leads to drastic simplifications. Accordingly, we will evaluate Eq. ( |35|) within a model in which 
each molecule has a narrow band of excited states. If we set Kij = + then the chiral terms which arise 

from performing the orientational average (see Appendix ^|) lead to the result 
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where now n is shorthand for rii,nj, D = D^-, D' = Dj/y, and the dots represent terms we dropped which do not 



contribute within the approximation we invoke in which the relevant excited states are strictly localized. (However, 
our results can be generalized to allow the excited states to extend over a small complex of atoms, if one simply lets 
the indices label electrons in complexes rather than those on atoms.) For localized excited states, all matrix elements 
are diagonal in their site indices. Nonlocal corrections to our results will be small in the parameter t/E, where t is a 
hopping matrix element which sets the scale of, the width of the band of excited states and E is a typical energy of 
the excited states relative to the ground state.L-3 However, it is important to check that these nonlocal corrections are 
not proportional to a lower power of l/Ru than the local ones we keep. An analysis of the relative importance of 
nonlocal terms is given in Appendix [d], where we show explicitly (albeit only for typical terms when both molecules 
are excited in the virtual state) that nonlocal contributions to the chiral interaction occur at the same order in I /Ru 
as do the local ones, but they are smaller by a factor of order t/E. This result justifies our subsequent neglect of 
nonlocal effects. 

If both molecules are in excited states in the virtual state "n," we may set i = i' and j — j' . If only one molecule, 
say the Jth one, is excited in the virtual state, then j and j' may be different. We will consider these two cases in the 
next two sections. 







IV. TWO-MOLECULE TERMS 



In this section we carry the sum in Eq. ( |38| ) over excited states \i,n; j, m) , in which molecule / is in state \i n ) with 
its ith atom excited to its nth state and molecule J is in state \j m ) with its jth atom excited to its mth state. Because 
we are dealing with localized states, these virtual states are obtained from the ground state only by interactions 
involving electronic charges on atoms i of molecule / and j of molecule J. Thus we no longer need consider here the 
presence of positive nuclear charges. Neglecting contributions of relative order (t/E) (as discussed in Appendix 
we may set i — i' and j = j' in Eq. (p8|), so that the contribution from virtual states in which both molecules are 
excited, indicated by the superscript (2), is 
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(39) 



We now evaluate this expression using the procedure of Eq. (t 
terms of Eq. (Efl) we write 



To illustrate the calculation for the first two 



Expanding the other matrix elements in accord with Eq. 
virtual state, we obtain 



:<D-% n m*j) D -\ • (4°) 

(pq) and recalling that both molecules are excited in the 



= (A^Az;)o„V; 2 ^ 3 + (Az>Az>)o n y'AV lz V> z W 3 } . (41a) 

Here we have dropped matrix elements like (Ar^ Q Ar^)o n which involve higher than dipole excitations and would 
therefore vanish for the p symmetry we have assumed for the low-lying excited states. In any event, since Ar ~ ao, 
the Bohr radius, this term would be smaller than those we keep by a factor of order (ao/R). To write the second line, 
one notes that D depends on and r'j only through z[ and Zj. Similar relations hold for the other terms in T\, e. g. 
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Thereby for the first two terms of Eq. (]39l) we obtain 
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Treating the other terms in Eq. ( |39|) similarly, we obtain 
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where p and v range over the labels x, y, and z of the local atomic excited p states and E vfl (i,j) is the energy of the 
virtual state relative to the ground state. (In principle, this energy can depend on the positions of the excited atoms. 
However, in our simplified treatment we will neglect such dependence.) In addition note that the expression given 

(2) i 1 

above for k\j must be averaged with respect to up and down directions of 7-th and J-th molecules, as in Eq. fll5p. If 
the excited states have a degeneracy with respect to spin, then the sum over \i and v should be extended to include 
a sum over spin indices. However, since singlet-triplet transitions are nearly forbidden, the multiplicity due to spin 
does not affect our results. Thus we obtain the result written in Eq. (|l6|). I— . .— . 

As discussed in the preceding section, our result is similar to that given by Van der Meer et al.jB and KatsJiil The 
important new aspect of Eq. ( fi"3| ) is that Kj] is expressed as a sum of contributions from pairs of atoms, one on each 
molecule. This formulation is consistent with the concept of local chiral centersH For L <C JZ, our expression for Pifj 
based on Eq. (fi"3|), when written in the form of Eq. (|l6|), reduces to that of van der MeercJ and KatsE-3 when My 
does not depend on i and j. However, when L is not much less than R, the fact that iSy involves an average over 
distances between atoms (rather than simply the distance between the centers of mass of the two molecules), leads 
to very different results. In any case, it is important to realize that My should be evaluated with respect to localized 
states, as is done here. 



A. Helical Molecule 



In this subsection we give a concrete evaluation of the above expression for two identical helical molecules. In the 
above formulae, position operators are given in the coordinate system fixed in the molecule while matrix elements 
are taken with respect to atomic p states which are referred to the principal axes locally defined for each atom of a 
molecule. Let us introduce the parametric representation of coordinates of an atom on a helical molecule: 



z = s, x = a cos(qs), y' = a sin(qs), 



(44) 



where q, the chiral wave vector of the helix, is defined so that a right-handed moleculeE3 has q positive. The locally 
defined principal axes for the zth atom at z' = s are chosen in the following way (see Fig. 0): 



e'4 = cos{qs)e' x + sin(qs)e^ 

e" y = c(- sm(qs)e' x + cos{qs)e' y - aqe' z ), 

e" z — c (—aqsm(qs)e' x + aq cos(qs)e' + e' z 



(45) 



where c 2 = [1 + (aq) 2 )] -1 . Here e" z is the tangent vector to the helix at z' — s, e" x is a unit vector along the radius of 
curvature at z' — s, and e'( y is the unit vector along the binormal or the third orthogonal direction.E2l We assume that 
the principal axes for excited p states coincide with these principal geometric directions. If we write e" M = Oi-^e'^, 



then the inverse transformation is e' = Oi^e" 



Note that the local axes are defined so that the matrix elements in Eq. 

(Mi|Ar^|0) = O i;p ,( M i|Ar;'|0) , 



t) are 



(46) 



where (/i,;|Ar"|0) is nonzero only if p = [i. Thus, in terms of local atomic coordinates we may evaluate Eq. ( |l8[ ) to 
obtain 



Mij = 3e 4 c 2 a 2 q 



x^Ax'nO) 2 (( Zj \Az>>\0) 2 /E xz - ( yj \Ay'!\0f/E xy ) 
c 2 [l 2a 2 q 2 }{y t \Ay>>\0) 2 ((z 3 \Az>>\0) 2 / E yz - (y^Ay'^O) 2 / E yy ) 
c 2 [2 a 2 q 2 ](z t \Az>m 2 ((y \Ay>;\0) 2 /E zy - ( Zj \ Az''\Q) 2 /E zz 



(47) 



This quantity can not depend on the locations of sites i and j because it is invariant against rotation about the long 
axis of the molecule and all locations on the helix are equivalent once end effects are neglected. Thus, neglecting end 
effects, we obtain the limiting results, 
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In both limits, the molecule is only weakly chiral, as we illustrate in Fig. ^. (To measure chiral strength the criterion 

of Ref. [2^ may be invoked.) 

(2) 

Now let us consider k\J as a function of the molecular length L. For simplicity we assume that the molecules are 
aligned exactly along their local nematic directions. Also we simplify the calculation by considering only the case 
when R/j is perpendicular to a/. Thus we will set 



a/-aj = l, (Dy • aj)(Dy • aj) = (zj- - 2j) , 



Then 
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(49) 
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(50) 



where N is the number of atoms in a molecule and p — N/L is the number of atoms per unit length in the molecule. 
For this simple calculation the average of Eq. (|l5|) is superfluous, so that nfj = k 2 J . The asymptotic result for 
L < R that 4j ~ 

can be seen in previous calculations. However, even in this limit, the fact that Kjj ls 
proportional to I? is not apparent from the previous results. To our knowledge, our result that k 2 ) ~ L/Rj for 
L > R is a new one. 

The macroscopic chiral wave vector Q and R are both taken perpendicular to the nematic direction. For QRu <C 1, 
we have a/ x aj-Rjj = —QR 2 cos 2 where <j>& is the angle between R/j and Q. Then the chiral energy per molecule 
from virtual states with two molecules excited, £^ 2 \ is given by 



\ YjPu) = -\Y,Q R2 cos2 0h4j 
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(51) 



In obtaining this result we approximated the sum over J by a sum over 7 nearest neighbors in the plane as specified in 
Eq. (|49]), so that cos 2 4>r — ► -|. ^From the discussion in Appendix |e] we are led to believe that the result of Eq. ( p3l| ) 
will not be seriously modified by taking a more realistic distribution of nearest neighboring molecules. We identify 
this result with the contribution to the torque field h in the Frank free energy from virtual states with two molecules 
excited: 
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where we took the volume per molecule, f2, to be fl — LR 2 for L > R and R 3 for L <C R- For L <C -R, M 2 ) <x (L 2 /R 9 ), 
consistent with the previous results of Refs. |l3| and |TJ, whereas for L > R, oc 1/R 7 . 

B. Numerical Estimate of the Macroscopic Pitch 

Now we want to estimate the value of the pitch using the result for £ obtained above. Intuitively one expects 
that the polarizability tensor will have its largest component tangent to the helix and that the anisotropy of the 
polarizability in the plane perpendicular to the tangent will be small. Essentially, we will attribute the anisotropy of 
the polarizability to the anisotropy in the excitation energy E a . Therefore, somewhat arbitrarily, we will take all the 
matrix elements like |(/i|Ar" Q |0)|, where p — x,y,z, to have same value, a a , where a a is of order the radius of an 
atom. We therefore parametrize the excitation energies in Eq. ( p7| ) as 

E x /E = 1 + \5 + T) , E y /E = 1 + ±<5 - 7} , E z /E = l-lS, (53) 

where E is the average excitation energy. Within our assumption of constant matrix elements the parameters S and r\ 
characterize the anisotropy of the excitation energy and through it the anisotropy of the atomic polarizability. When 
this anisotropy is small, we find that 



where 



and 



^ = (1 + aV) (55) 



G(5,r,,aq)=( a ? {6 2 j ) *(ag) ■ i5(i) 



1 + a 2 q 2 



(2) 

The corresponding results for Kjj are 



K?} = - 3e4 °l a £ L2 G(6,r,,aq), L«R (57a) 
45*Ata?L L ^ R 



128ER 7 

Thus M is quadratic in the anisotropy of the polarizability and 

h(2) = _ 3^aiap^ G{§T]aqh L<<R (5ga) 

W - 45 Sy ^"^' L> R . (58b) 

This conclusion is a natural one: surely the torque field must disappear when the anisotropy of the polarizability is 
turned off. Also, when E z = E y (i. e. when 5 = rj), the chiral constant a 1 vanishes. To see that note that when 
E z = E yi one of the principal axes for each atom can be taken to be perpendicular to the axis of the helix, in which 
case the matrix elements appearing in a' are invariant with respect to the mirror operation z' — > — z'. To illustrate 
the dependence of hf^ on the molecular chiral wavevector q, we show in Fig. || G(5,i],aq) versus aq for fixed values 
of 5 and rj. There one sees that is maximal for aq of order unity and decreases rapidly away from this maximum. 
Of course, an experimental test of this dependence is difficult since varying q at constant p involves structural changes 
in a molecule. To treat small chirality we take aq = 1/3 (or aq — 3) and we set a a — 1A, E = 8 eV (these parameters 
correspond to an atomic polarizability a — 2e 2 a 2 /E = 27a,g), a = 7. 5 A, 7 = 6, L = 200 A, R = 20 A, p — 3 A , 
S = 1/5, and r\ = 0. With the volume per molecule, ~ LR 2 , the chosen values of the parameters correspond to 
volumetric density of molecules of about 40% and a dielectric constant, e = 1 + 4irapL/£l w 1.3. Then the torque 
field is approximately h — 4.5 x 10~ 4 (dyne/cm). If now one takes the Frank constant K2 to be 10~ 7 dyne, then the 
macroscopic pitch of the liquid crystal will be P — 2-k/Q — 2ttK~2/ h = —14// (or 28p for aq — 3). If we had taken 
S = 3/10 and ij = 0, then the pitch would be — 4.5/t (or 9.0/x for aq = 3). 
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It may be seen that the computed pitch is longer then one usually finds experimentally for a system consisting of 
molecules of the above size. There are two possible explanations for this discrepancy. First of all, our approximations, 
although improved over previous ones, may still not be sufficiently accurate. For example, for two helices of radius 7.5 
A at a center-to-center separation of 20 A have their nearest groups separated by only 5 A. Under these conditions, 
the expansion in terms of even the transverse coordinates of the atoms may not be rapidly convergent. The second 
possible reason for the discrepancy between calculated and observed pitches would be that an explanation of the pitch 
of cholesterics requires consideration of steric interactions. We are presently considering how our arguments might be 
improved to discriminate between these two explanations. 

If one can find molecules for which quantum chiral interactions considered in this section are dominant, then the 
following remarks are relevant. Notice that for helical molecules the torque field, h, can have either sign in both 
the large q and small q limit, depending on the. signs of (5 — rf) and (5 — 377). This is in contrast to the situation 
for steric interactions, for which it is believed" that the contribution to h from the repulsive (i. e. steric) chiral 
interaction between molecules is negative for small q and is positive for large q. Helical molecules which do not follow 
the sign prediction for h due to repulsive steric interactions might constitute examples of molecules for which the 
quantum dispersion forces dominate the chiral interactions. In general, the density dependence of the quantum and 
steric contributions to h will be different. Thus, if these two mechanisms compete, it is likely that the sign of h could 
depend on the density. 



V. ONE-MOLECULE TERMS 



In the model of a molecule considered before we supposed it to consist of He-like atoms. In reality one would expect 
the outer electronic shell of atoms to be deformed by the interaction with nearest neighbors. In general, constituent 
atoms or complexes will possess a dipole moment. Hence it is of interest to consider the situation when one of the 
molecules is in its ground state in the virtual state of two-molecule system. Up to now this case was ignored, although, 
as we shall see, it may play a significant, if not dominant role. 

^From Eq. (^) we obtain the following expression for the additional contribution, denoted k\ j , to Kjj from itirtual 
states in which only one molecule is excited. We still invoke the approximation of localized excited states.EJ But 
then terms in which only molecule J is excited require evaluation of £ij- t i>j> n with j = j' , but i and i' are arbitrary 
and similarly when only molecule I is excited. For a molecule in the excited state we use the same approximation 
as before, again expanding the denominator with respect to Ar to get a nonzero matrix element. For the molecule 
which remains in its ground state in the virtual state, one has to include both signs of charge at each site. Thus (see 
Appendix |]) we find that 

r 2 1—1-1 —r ->\ (Pi'j ' a J) 

«>j =6 > e q l q v [x l y i , - x v yA _ 3 _ 5 — 

/ — n n 

i,i'£l;j£j ij i'j 

x ^^(ir^^-IA^IO) 2 - |< Mj |A^.|0> 2 - i< Mj lA^10) 2 ) ■ (59) 

In Eq. (|59j ) we sum i and i' over all the charges in a given atom, in which case qix^ is replaced by p' xi , where 
is now the expectation value of the dipole moment of the ith atom, in its ground state, so that i and i' from now on 
refer to atoms, whereas j will still label electronic charges. Then the preceding equation can be reduced to 

IJ — / ylPixPi'v Pi'xPiyl — 3 — 5 
ii'j ^ij i'j 

x e^M^m^W - !<M;|A4|0> 2 - i( Mj -|A^|0) 2 ] (60a) 

= WiW 2 , (60b) 

where W\ is the factor on the first line of this equation and W2 that on the second line. In writing this result we 
assumed that for typical atoms, i, one has Pi X ~Zi' S> Piz%i'- Once again, in this expression one has to carry out 
averaging with respect to independent up and down orientations of both molecules. But this average turns out to be 
superfluous for the model of a helical molecule which was introduced above. 

As in Eq. (f45|), we introduce components of the atomic dipole moment with respect to the principal axes of the 
atom, in which case we have 
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Px = Px cos(gs) — cp'y sin(qs) — caqp" sin(gs) 
v' y = Px s'm(qs) + cp'y cos(qs) + caqp" cos(gs) 
p'z = -caqp'y + cp" z . 



(61) 



The component p" x is essentially the radial component pf the atomic dipole moment and is nonzero for helical geometry. 
For instance, for the molecule TMV, shown in Fig. pp3 this radial component may be appreciable. In such a case we 
write 



2 / / 

Pii' — PixP 



txt-vy Pi'xPiy = [Px + c (p'y + aqp") } sm[q( Sl/ - Si )} 
= plsm[q(si> - s^] . 



We now substitute this form into Eq. (pOal) and assume perfect alignment as in Eq. 
symmetrized and we write W\ = p\X\L j R? , with 



X 1 (q,L) = 3j2 



1 - L 2 (s~ 2 + S.S,, - Sjjsj + §i,)) 
[1 + L 2 {s J - .? 4 ) 2 ] 5/2 [l +L 2 (Sj - S e ) 2 } 5 / 2 



(Si 



I sin [<?(,' 



where s = s/L, q = qL, and L = L/R. To evaluate W2 we again invoke the model of Eq. 
small anisotropy, Eqs. ( p6|) and (Q) enable us to write 

W 2 ee e^i^K^IA^lO) 2 - i<^|A.T;|) 2 - i(^|A y ;|) 2 ] 

A" 
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(62a) 
(62b) 

Then the summand is 

50] . (63) 
(HH), in which case, for 
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where ^(aq) is defined in Eq. (|55|). Using the asymptotic evaluations in Appendix |g|, we thus have the results 

4j = a ERS V(aq)<t>(±qL) , a « L « R ; 



ie 4 a 2 d 2 p 3 qL 
ER 1 



V(aq)If(qR) , L » R ; 



where (i is the effective size of the dipole moment: po = e -d, 

<f>(x) = -{3/2)(d/dx)[(sinx)/x] 2 

and 



I n (qR) = 



e 2 



(64) 

(65a) 
(65b) 

(66) 

(67) 



Now we evaluate h following the procedure of Eq. ((5^) in terms of the chiral energy per molecule £^ due to 
one-molecule effects: 
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Using the evaluations of Appendix |g], we obtain the asymptotic results, 

7e 4 a 2 d 2 LV 3 



^(aq)X 1 (q,L) . 



4ER 9 
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ER 4 



V(aq)<t>(hqL) , a<i<i?, 



V(aq)I*(qR) , 



(68) 

(69a) 
(69b) 



Here again we see from the appearance of ^(aq) that chirality requires a nonzero anisotropy of the polarizability 
characterized by S and 77. Since the factor ^(aq) also appears in Eq. (El), we see that the critical value (if any) where 
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h changes sign as q is varied, is only determined by the geometry at least within our simple model. For concentrated 
systems, the limit L ^> R is the most relevant and for this case Fig. fj] shows how depends on the molecular 
chirality q, when the length of the molecules and the density of atoms p are fixed. Note that the variation of 
with the molecular chirality q strongly depends on details of molecular geometry since only a fixed number of atoms 
is allowed on a helical thread. Figure [?] shows that for aq of order or less than unity, where is appreciable and. 
may give a short pitch, ft/ 1 -* is positive, whereas for steric interactions h is believed to be negative for small gEta 
Since increasing the density probably causes steric interactions to dominate, it is possible .that the combination of 
these two mechanisms could cause h to change sign as a function of density or temperature.Ej In Fig. || we show the 
behavior of the quantity Y\ = ^(aq)Xi(q, L) as a function of L for R = 20A, p — SA^ 1 , S — 1/5, and rj = for two 
fixed values of the molecular chiral wavevector, q. In particular, it is noteworthy that for large L, Y\ (and therefore 
h^) is independent of L. To get some idea of the relative importance of and h/- 2 \ consider their ratio: 

L < R ; (70a) 
L > R . (70b) 

One sees that even with d/a a as small as 0.03, this ratio can easily be of order unity. 

To numerically estimate the pitch arising from the considered interaction we will take parameters of a system 
and constituent molecules chosen in preceding section. Then, if aq = 1/3 one finds = 0.5(d/a ) 2 dyne/cm. 
If molecules posses a local dipole moment, the resulting dipolar interactions may lead to strong biaxial correlations 
between neighboring molecules. Using the evaluation of the dipolar interaction energy in terms of the integral analyzed 
in Appendix |G], we estimate that the order of magnitude of the dipole-dipole interaction to be Vdd ~ (d/ao) 2 ■ 10 5 K. 
So if we suppose that the biaxial correlations due to dipole-dipole interaction among molecules is negligible when it is 
less than 100K then one must have (d/ao) 2 < 10~ 3 . At the upper limit of validity of our calculations (d/ao) 2 = 10~ 3 
and the macroscopic pitch due to will be PW — 2ttK2/ 'h^ — 12.5 p. As the density of local dipoles is increased, 
the macroscopic pitch becomes smaller. For instance, if we set d/ao = 1/3, we get a pitch of order O.lp, although 
this estimate will be significantly modified by biaxial correlations, which have been neglected in our treatment. Since 
/jW ^ p 3 and the dipole-dipole interaction is proportional to p 2 , it is conceivable that for much larger molecules h^> 
could be significant without the dipoles being large enough to induce long range biaxial order. Finally, when aq is 
larger than unity (as for TMV), this mechanism leads to a very large pitch for almost any choice of parameters. As 
mentioned in Sec. IVB, it is possible that larger values of the pitch would be obtained if the role of the transverse 
were treated exactly rather than by an expansion. 

VI. CONCLUSION 

Here we put our work into the context of current research and record our conclusions. 

1) We introduced a simple model of localized polar excited states that enabled us to make an explicit calculation of 
the chiral interaction, kijsli x aj • Rjj, between molecules / and J due to quantum charge fluctuations analogous 
to those responsible for the R~ 6 dispersion interaction between neutral atoms. We identified two distinct physical 
effects depending on whether one or both molecules were excited in the virtual state of the two-molecule system. In 
implementing this calculation we used a modified multipole expansion in which only coordinates transverse to the 
long axis of the molecule were expansion parameters, so that we could treat long molecules which usually are the 

(2) 

building blocks of liquid crystals. The contribution, js^j ' , to K/jriom virtual states with both molecules excited has 
a form similar to that found by van der Meer et alii 3 ! and Katslij. For a helical molecule of length L we find that 
nfj oc L 2 /R s for L <^ R and Kj) oc L/R 7 , for L > Rij. The contribution, nf] to Ku from virtual states with only 
one molecule excited is usually only dominant when the local atomic dipole moments are large enough to give rise 
to significant (possibly long-range) biaxial correlations. Both mechanisms give rise to a chiral interaction between 
a chiral molecule and an achiral one that has a local anisotropic polarizability. Our formulation leads to numerical 
estimates of the pitch which are larger than that found in many cholesterics. Whether this discrepancy is an artifact 
of the expansion in transverse coordinates along with a disregard of biaxial correlations between molecules or is an 
indication that steric rather than quantum interactions are the microscopic origin of macroscopic chirality is not clear 
at present. The role of biaxial correlations between molecules will be considered elsewhere.Ej 

2) We evaluated kjj and the torque field h for helical molecules as a function of the wave vector q which describes the 



fcW (2L 2 d 2 p\ fl + a 2 q 2 



M 2 ) \ 3a 2 a J \ aq 
( 1024R 3 d 2 p 



yS-nr^iiqL), 



V 45irala 2 



(l + a 2 q 2 )(5- V )-'lf(qR) 



1G 



chiral structure of an individual molecule. We found that the sign of h depends on the details of the anisotropy of the 
local atomic polarizability. For instance, for (aq) 2 <C 1, the sign of ftA 1 ' (the contribution to h from virtual states in 
which only one molecule is excited) has the same sign as S, the local anisotropy of the polarizability. One expects S to 
be positive because presumably the polarizability along the tangent of the helix is larger than that along perpendicular 
directions. This sign of hP^ is opposite to that expected from steric interactions.^ As for steric interactions, one 
expects to change sign as q is increased, but our calculations indicate that this only happens when is so 
small that it is hardly likely to be the dominant mechanism for macroscopic chirality. When 6 is positive and large, 
the sign of the two-molecule contributions to h is negative for small aq and positive for large aq, just as expected for 
steric interactions. However, our calculations indicate that normally h^ 2 ' is not significant. 



3) Here we calculated the effective chiral interactions by averaging the orientation of the molecule over configurations 
with the long axis fixed. Even within mean field theory, wherein each molecule is described by a single-molecule 
orientational distribution function of the three Euler angles, the only required symmetry in the locally nematic state 
is that it be invariant against rotations about the nematic axis. As discussed in Appendix B, this requirement still 
permits biaxial contributions to the orientational probability distribution which we neglected. 



4) These calculations suggest some general observations. First of all, the interaction from virtual states with two 
molecules excited, give rise to a two-point chiral interaction in the form of an integral over the long axis of each 
molecule. This result gives a formal justification for introduction of a chiral interaction between "chiral centers" on 
one molecule with a center of local anisotropic polarizability on another molecule. However, this same characterization 
does not apply to the mechanism involving local permanent atomic dipole moments. The dipolar mechanism leads to 
an intrinsically three-point chiral interaction of a type which, as far as we know, has not yet been proposed, ft would 
be interesting to observe such an interaction for helical molecules which have a local radial dipole moment. 



5) Our calculations can potentially be generalized in several directions. For instance, there seems to be no reason why 
our results can not be taken over immediately to discuss the interaction between flexible polymers. There the average 
over spinning (within a tube surrounding the convoluted polymer shape) can still be taken. Then in Eq. (16) one 
would replace a/ by its local value at atom i. Our calculations can also be applied to liquid crystal systems containing 
a mixture of chiral and achiral molecules. There one has two types of interactions to consider. The first of these is 
the interaction between adjacent chiral and achiral molecules to which the results of this paper apply directly. The 
second is the interaction between more widely separated pairs of chiral molecules. For this interaction, our result for 
kjj ought to be multiplied by e~ 2 , where e is the static dielectric constant. 
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APPENDIX A: QUANTUM AND CLASSICAL AVERAGING 



The energy of interaction of molecules I and J averaged over their rotational motion when expressed in terms of a 
multipole expansion is of the form 

[Ulj] av = ^ Aa 1 ,a 2 ,...a n ;/3 1 ,l3 2 ,...l3 m (I, J)[M ait a 2 ,...a n (pi)]a,v[N0 u 2t ...0 m (pj)] a . v , (Al) 
M{a},N,{/3} 

where [ ] av indicates an average over orientations and M and N are tensors of arbitrary rank which are functionals 
of a density on the molecule in question. For classical two-body interactions these tensors are multipole moments of 
the form 



. (pi) = J dvp!(Y)r ai r a2 . . .r an . (A2) 



For classical two-body interactions these tensors are thus linear functions of the density, so that the orientational 
average of the tensor is the same as the tensor evaluated for the orientationally averaged density: 

[M a (p)]av = M ai (Mav) • (A3) 

This means that classically the interaction averaged over the orientational motion of molecule /, say, is the same as the 
interaction would be for a molecule having the average (over orientations) shape. Thus, classically, spinning a chiral 
molecule leads to two-body interactions characteristic of a uniaxial, i. e. achiral molecule. Quantum mechanically, 
the situation is different, because in second order perturbation theory the tensor M, say, in Eq. (Al) is a bilinear 
function of the density p(I) of the form 

M ai>a2 ,... an = J p(r)dr J dr'p(r')T(r, r')r ai . . . r ak r' ah+i . . .r' an , (A4) 

where T(r, r') depends on the spatial correlations of the important excited states, and Eq. ( |A3| ) is incorrect. In other 
words, the nonlinear fluctuation of the electric field of a molecule due to quantum fluctuations has a chiral component 
that survives an average over rotations and thereby distinguishes between right-handed and left-handed molecules. 



APPENDIX B: BIAXIAL ORIENTATIONAL CORRELATIONS 



If the Euler angles are taken to represent the orientation of the molecule with reference to axes fixed in space such 
that the z-axis coincides with the axis of nematic order, then the probability distribution for the orientation of a 
single molecule must be independent of a. If the distribution is also independent of 7, then it means that for each 
value of P, the molecule spins with equal probability through all angles about its long axis. However, if we have 
correlations between j3 and 7, we can have a distribution like that describing the orientation of the moon in which 
7-/3 assumes a fixed value. For a molecule, this distribution is depicted in Fig. Uu. 



APPENDIX C: ORIENTATIONAL AVERAGES 



In this appendix we evaluate the orientational averages (indicated by brackets, [ ] av ) of the terms in Eq. (B50. In 
this calculation, we should keep in mind that we only need keep terms which include one antisymmetric tensor. Also 
only averages of even numbers of powers of components of p i are nonzero. Finally, terms obtained by interchanging 
the indices i and j (labeling atoms on different molecules) can be included implicitly. With these understandings we 
use Eq. (|l^) to write 



T 1 = [(p ij -r>D- 3 ) n0 (p ij -T>D- 



)0n\av 



— 2[(pi a D a D )no(PipDpD )on]av 

= ea/3 7 a/ 7 e MI , z (r. D a D~ z ) nQ (r' iv D D~ 



)0n 



[R x ai • aj] 



(r' ifi ZjD- 3 ) nO (r' iu D- 3 ) 0n - (r' ilx D-*) n0 (r' ivZj D-*) Qr , 



/ /l 

2 [R x ai • aj] (r'^D-^^D- 3 )^^ 



(CI) 
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In terms involving four powers of transverse components, contributions at the order in Rjj which we need require 
that two components refer to atom i and two to atom j. Thus 

T2 = -| [([ Pij ■ D] 2 ^- 5 )„ (p^- 3 )o„] av 

= -3 [([ Pl ■ D] 2 C 5 ) n0 (p^-3 )on ] av _ 6 [ ([p . . D][p . . D])„ ([p j ■ p,]^- 3 ) „] av . (C2) 

The first term gives rise to no antisymmetric terms and can be dropped. The second term leads to 

T 2 = -3 [(r^D aP]0 D p D- 5 ) nO (r^p ja D^) Qn ] av 

+ 3 [( r i f * D aai a PjpD l 3D~ 5 ) n o(r' ifl a Il p n D~ 3 ) 0n 

-3 e^ z e al3p ai p ((r^D a p J f l DpD' 5 ) n0 {r' w p : i J D~ 3 ) 0n ] iiv . (C3) 

The first term gives zero antisymmetric contribution. The second and third terms give identical contributions. So 

T 2 = 3 [R x ai • a, 7 ] (r^r^[D ■ ajp-^o^r^-^^^ . (C4) 

Likewise, keeping only relevant terms, we write 

T 3 = | [([ Pij ■ B} 2 D-% lQ ([p l3 ■ T>] 2 D- 5 ) Qn ] av 
= [|([ Pl .D] 2 ^ 5 )„ ([p r D] 2 ^ 5 ) „ 

+9([ Pi ■ D][ Pj • D]£- 5 )„o([ft • mpj ■ D]^V] av . (C5) 

The first term leads to zero antisymmetric contribution. In the second term there are two equal contributions, one 
from taking the antisymmetric term in the average over p iy the other from the antisymmetric term in the average 
over Pj. So we write 

T 3 = 9e^ z [{r^PjpD a DpD- 5 ) n0 (r / w p j sD 7 D s D- rj ) 0n e aip a Ip ] a:v 

= ^e^ z (r[ p r[ JT D a DpD^) n0 (r[y jT D 1 D s D^ 5 )( 3n e aip a I p{5ps - ajpa JS ) . (C6) 
Then, using the symmetry between the two matrix elements, we have 

T 3 = 9 [R x ai ■ blj] (r' ip r' ju D a D^ 5 ) n o(r' lll r' jT D 7 Dsz' i D^ 5 )onei'Tz(S ai - a/ Q a/ 7 ) ■ (C7) 
We set z[ = "z^ + Az^. The term in vanishes. Thus 

r 3 = 9R x ai • 2Lj¥ ip {r[ pi r l ju D a D~ 5 ) n n(r' ]T D 1 D8{&.z l i )D~ 5 ) 0n e VTZ {8 ai - a/ Q a/ 7 ) . (C8) 

The matrix elements are symmetric functions of p and r. So the antisymmetry of the e tensor causes this term to 
vanish. At higher order in Rjj there would be nonzero contributions from this term. But at the order we consider 
there are none. 

The remaining terms in Eq. ( |35| ) vanish for reasons similar to those which made T3 vanish. So the only contributions 
that survive are those written in Eq. (|38|). 

APPENDIX D: NONLOCAL EFFECTS 

In this appendix we discuss nonlocal corrections contained in Eq. ( |38| ) from terms where i 7^ i' and/or j ^ j'. 
Rather than give a general argument, we will illustrate the nature of the argument by considering specifically the 
nonlocal corrections to the first term in Eq. (|38|). For this purpose we assume that the "unperturbed" energies E riinj 
can be obtained from a Hamiltonian of the form 

H = Ho + Vhop 

^0+ E Wo'fl > (di) 

where TIq is completely local: 
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Ho = Yl \ ia ) E i 



i a 



(D2) 



We assume the states to be strongly localized so that \t^\ -C Ef for all indices. 

Now we consider the contribution, Tb> to kjj — \ + Hji] from the first line of Eq. (Bq). Thus we write 



T = -e 4 



\n a D^-r llp z r D/ r \ 



(D3) 



In this appendix all coordinates are taken relative to axes fixed in the molecule. Thus, r^ 1 here denotes what we called 
(r')^ in the notation of Eq. (p~T|) . For simplicity we consider here only the contribution from virtual states in which 
both molecules are excited. In that case, the sums are only over electrons. 

Now we expand the matrix elements according to Eq. (|36|), as was done in Eq. (pll|). Thereby we get the 
corresponding contribution 5kjj as 

Shij = (0|Ar^Ar^-Ari' CT Arj T |0) 



\Ar i( .Ar jri -Ar v<J Ar jT \0) 



f 3 _3 3 



(D4) 



Here we dropped terms of order 1/R^j. In evaluating the gradients, note that Dij depends on (tj) only via Zj, (zj). 

3 

Thus V^V jr/Dij is only nonzero for £ = r/ = z. 

The terms of greatest interest are those of order l/Rjj, because such terms are of potentially lower order than the 
local terms we kept of order 1/R^j. These leading order terms are 

Sku = {Q\Ar il± Ar jz -A n ,„Ar jz \0) (d;,],V 3Z D^- -D^V r Z D^\ e „„ a 



8 

{Q\Ar itJ ,Arj z -Ar Vv Ar jz \Q) 



f(Zi, Zj, , Zj' ) f(Zi' i Zj' , Zi, Zj~) 



(D5) 



where / ~ l/Rjj. Note that when the states are locali zed, i. e. when i = i' and j = j', the factor in large square 
brackets vanishes. Now consider expanding £ as in Eq. (Dl), so that 



-V h 



C TP 'XJ ' TP "\J ' Il0 P 77" 0_/ ' JP <\J ' TP *\J ' TP "\J 

c h,o — no &o — no j^o — no &o — no &o — /TO H/Q— no 



lop 



lop 



(D6) 



where Eq — TLo ~ E, where E, the typical excitation energy, is much larger than i, a typical hopping matrix element. 
This equation implies that when it requires m hops for an electron to move from site i to site i' and n hops for an 
electron to move from site j to site j' , then the matrix element will be of relative order (t/E)( m+n \ Thus 

A/ = f{ Zi ,z h z v ,zy) ~ f{z^z r ,z uZj ) ~ (t/E)Vf ~ (t/E)( Zi - z i+l )/R s u . (D7) 

We see that the ratio of this nonlocal contribution to the local contribution of Eq. ( |38| ) is of order A/ /(r±/ Rjj), 
where r± is a typical value of Xi or y,^ This ratio is thus of order (t/E)(zi+i — Zi)/r±. Normally (z^+i — Zi)/r± is of 
order unity, so indeed the nonlocal contributions are of relative order t/E and can reasonably be neglected. 
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APPENDIX E: POSITIONAL CORRELATIONS 



In this appendix we consider how energy of interaction for the system of molecules is effected by the relative 
distribution of molecules. A simple way to address this issue is to evaluate the chiral interaction as a function of 
Zu = Z, the z-component of R/j. We assume that it suffices to do this for helical molecules, in which case the 
calculations can be done explicitly. Previously we had set Z = and had considered the contribution to the torque field 
from a shell of six neighbors taken to lie in the equatorial plane. Here we show numerically that this approximation 
is quite reasonable. We study the dependence of Su = Y]jj Sij on Zjj. We still assume perfect nematic order, so 
that aj = a. j = e z . Then the sum in Eq. ( pp| ) becomes 

Su(z) = J2lR 2 + (z + Zl - Zj ) 2 ]- rj [R 2 ~(Z + Zi- Zj f] , (El) 

ij 

and we see that 

_ Su(Z) 

4 2 j(o) ^(o) ' [tjZ) 

For the one-molecule terms we similarly note that the ^-dependence in Eq. ( p9[ ) is reproduced by writing 

— n n 

ii' j ij i'j 



E 



p 2 (Z + 4 - z' 3 ) sin[<j(4 - z'i)] ^ 



- [K 2 + {Z + ^-z') 2 ] 3 / 2 [i? 2 + (Z + z' t , -z') 2 ] 5 / 2 



in the notation of Eqs. ( |59| ) and ( |6 

These results allow us to compute the ratio k^J (Z)/k^J (0) which is shown in Fig. ||for n = 1 and n = 2. This 
result is representative of the situation for a wide range of parameters. As one might expect, the contribution to the 
torque field decreases strongly as \Z\/ L increases towards unity. Accordingly, the approximation of including only the 
effect of equatorial neighbors is a good one. 



APPENDIX F: CONTRIBUTIONS TO H (1) 



In this appendix we discuss the evaluation of the one-molecule contributions to kjj. We consider the terms in the 
last two lines of Eq. (|38|). We will analyze the one-molecule contributions which arise when i — i' but j and j' are in 
general different. In the intermediate excited state only atom i is in an excited state. Atoms j and j' remain in their 
ground states. Calling this term T we write 

T = -3 g <-V/, 0.//V /),/!,•;, 'D, ; • a, />,/'(> + . . . , (Fl) 

i,j,j' 

where the dots denote the three additional terms required to make the expression be rotationally invariant. (These 
can be reconstructed at the end of the calculation.) Using the expansion of Eq. ([56]) we have 

T = -3e 2 ^ mi , (0| [Ax^r.f + x[&z[{V lz D-?)]y' r ~ 

ijj' 

x [Aa-JD^ • ai Dr/ +x' i Azy iz [CD lJ ■ a z )£>r. 5 ]] 2j.|Q) + . . . 
w -3e 2 ^ my (OIA^^A^IO)^^.,^ (D^ • a,)^ 5 + . . . . (F2) 

ijj' 

Now we carry the sum over j (j 1 ) over the charges that comprise the dipole moment pj (p,-/) on atom j (j') to get 
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r = -3e 2 ^(0|A^iA^|0)[p^^+^,(V^^?)] 



'jj' 



x [p' xj (Dij ■ ai)D y B +p' zj W' j ,V jz [(D ij ■ m)D7. 5 ] 
-3 e 2 ^(0|Ax^A^|0) P ; / p^^f(D y - • ai )D^ 5 



ijj' 



where now j and j' refer to atoms. Restoring the additional terms to preserve rotational invariance we obtain 



T = -3e 2 ^[<0|A^A^|0) + (OIA^A^O)]^,^ - p' X] ,p' yj \D^ (T> 13 . ai )D: 

ij'j' 



(F3) 



(F4) 



When the indices are relabeled, this result reproduces part of Eq. (pit). 



APPENDIX G: EVALUATION OF INTEGRALS IN SEC. V 

In this appendix we evaluate the integral X\ in Eq. ( |63| ) and an integral needed to evaluate the dipolar interaction 
energy between two long helices. 



Consider the asymptotic evaluation of Eq. (63), firstly, in the limit L 3> R- End effects can be shown to be 
negligible, in which case the final summation (over Sj) introduces a factor of N and one sets Sj = 0. Also we consider 
only the continuum limit in which the sums are replaced by integrals. One can show that correct to leading order in 
L , the limits on the integrals can be extended to ±00. Thus we have the asymptotic result 



X!(g,L)~3iV £ 



ds I ds' 

-00 



I - L 2 ss' 



(I +L 2 S 2 ) 5 / 2 (l+ZV 2 ) 5 /2 

For each of the two factors in the denominator we introduce the representation 

1 



(s' — s) sm[qL{s' — s)] 



P 



-5/2 _ 



x 4 e"2^ dx 



(3V2^) J- c 

Then the integrations over s and s' can be done analytically and eventually one finds that 

X 1 (q,L - 00) = 8{pR) 3 (qR)li(qR) , 



where I\ is defined in Eq. ( |67| ) of the text. 
The limit L — > is trivial. We find that 



X 1 (q,L^Q) = N 3 ( f>(±qL) , 



(Gl) 



(G2) 



(G3) 



(G4) 



where </>(x) = — (3/2)(d/dx)[(sinx)/x] 2 . 

Finally we evaluate the dipolar interaction energy E dd between two long helical molecules, a and b, separated by a 
distance R along the x-axis. We assume that the radius of the helix is much less than R. In this limit, in terms of 
the atomic dipole moments we write 



E dd = p 1 



L/2 



rL/2 








/ dz b 


R 2 + z 2 ab 


-3/2 | 


p"cos(gz Q + 4> a ) 


/-L/2 









dz a 

-L/2 

p" x cos(qz b + (p b ) - cp'y sin(gz b + 4> b ) 
p" x sin(qz a + 4> a ) + cp'y cos(qz a + (f> a ) 
p" sin(qz b + 4> b ) + cp'y cos(qz b + 4> b ) 



3R 2 



R 2 + *l b 



(G5) 
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where py 



Py + a QPzi 0a (</>&) is the angle of rotation of molecule a (b) about its long axis, and z a b = z a — Zb- Here 



we did not include terms involving p' z which either are independent of both angles <p a and <p b or vanish in the limit 
L — > oo. In that limit we only need to keep terms which depend on z a b, in which case we have 



Edd 



1 2 2 t 
2P P0 L 



COs(qz ab + (j) ab ) 



2z 2 



R 2 



{B? + zl b )W 



where cf> a b = 4>a — <fib- Using Eq. (G2) we obtain the final result 

E d d = -Lp 2 plq 2 I_i cos0 q6 = -^V d dCOS(j) a b ■ 



(G6) 



(G7) 
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(a) (b) 



FIG. 1. Left: Molecule-fixed coordinate system, defined by the unit vectors e^. Right: Definition of the Euler angles a, (3, 
and 7 which take the space-fixed axes e^, e y , and e z into the molecule-fixed axes, e' x , e' y , and e' z . Note that a and j3 are the 
usual spherical angles which specify the orientation of the long axis of the molecule, e' z , with respect to the space-fixed axes. 
The third Euler angle 7, not shown here, is the angle of rotation about the z' axis which brings the x and y axes in coincidence 
with the x and y axes fixed in the body (respectively e' x and e^). 




FIG. 2. Space-fixed coordinate system, showing the displacement, Hi, of the 7th molecule and the displacement, r; of the 
ith charge of the I molecule relative to the center of the molecule. 
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z 



FIG. 3. Local atomic coordinate system, denned by the unit vectors e^, showing that the local excited p states define the 
orientation of the local axes. Here e" is the unit vector tangent to the helix, the unit normal, e", lies along the radius of 
curvature, and the binormal unit vector e'y is the third member of the triad of mutually perpendicular unit vectors. 





a e" 



FIG. 4. Locally defined principal axes for weakly chiral molecules with large q (left) and small q (right). No te th at the axi s 
nearly collinear with the long axis of the molecule is the j/-axis for large q and the z-axis for small q. In Eqs. (48a) and (48h) 
the anisotropy of the polarizability needed is with respect to the long axis of the molecule. 
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0.010 




0.006 



FIG. 5. The function 



G(5, rj, aq) = ^ i+'Sq' 1 J — r l)^'{ ac l) versus aq for 5 = ~ and 77 = 0. 




FIG. 6. TMV, taken from Ref. |35| We indicate a possible axis along which the dipole moment of each complex might be 
oriented. In the situation shown here, the largest component of the dipole moment of the complex is radial. 
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FIG. 7. Yi(qL,L/R) = ^f(aq)Xi(qL,L/R), with *(ag) and Xi(q, L) denned in Eqs. (||) and (||) respectively, versus qL 
for L — 200A, R = 20A, p = 3A _1 , 5 = 1/5, and 77 = . According to Eq. ( |68| ) the quantity plotted gives the dependence of 
the torque field on the chiral wavevector of a molecule q. Note that the molecule is achiral if either q — > or q — > 00. 
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FIG. 8. The factor Yi = 1 i/(aq)Xi(qL, L/R) as a function of L for two values of the molecule wave number q — 0.0444A 1 
(plot 1) and q = 0.1333A -1 (plot 2) with R = 20A, a = 7.5A, S = 0.2, and r? = 0. According to Eq. @ when L > R (so that 
f2 = LR 2 ) the quantity plotted gives the dependence of the torque field /2/ 1 ' on the length L of a molecule. 
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FIG. 9. The ratio KjJ (Z)/k^j (0), where Z is the z-component of Hij for molecules of length L = 200A and intermolecular 



separation R = 20 A. For n = 1 we show essentially indistinguishable curves for q ■■ 
n — 2 this ratio does not depend on q. 



0.0444(A)" 1 and for q = 0.1333(A) -1 . For 




FIG. 10. Distribution for which 7 — j3 is fixed. 
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